Vortex, Defect, and Monopole Gauge Fields for Nambu-Goldstone Systems 
and their Phase Transitions 



or 

Theory of Fluctuating Multivalued Fields 



H. Kleinert 

Institut fur Theoretische Physik 
Freie Universitat Berlin 
Arnimallee 14 D - 14195 Berlin 



Abstract. We survey the various uses of singular gauge fields which render 
an ideal description of ensembles of line- and surface-like excitations and 
explain phase transitions with the help of simple quadratic actions. Near 
a transition, the singular gauge fields can be transformed into ordinary 
complex fields called disorder fields which differ from Landau's order fields 
by developing a nonzero expectation value in the high-temperature phase. 



1. Introduction 

Many physical systems possess massless excitations. These dominate their statistical 
properties at low temperatures. The specific heat, for example has the characteristic 
temperature behavior 

C~T D (1) 

in D space dimensions (see Fig. 1 for solids and Fig. 2 for superfluid 4 He). 

The origin of massless excitations can often be explained by a spontaneously 
continuous symmetry of the action. Then the excitations are called Nambu-Goldstone 
(NG) modes. Examples are sound waves in superfluids and solids and magnons in 
ferromagnets. 
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□ Ag 215 K 5.61 K 
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Figure 1. The specific heat of various solids. By plotting the data against the 
ratio TV On. where On is the Debve teniDerature, the data fall on a universal 

temperature T m . 



free Bose gas 




Figure 2. The specific heat of superfluid 4 He. The peak is caused by the 
superfluid-normal transition. A T 3 -behavior is clearly visible under the peak. 



In superfluid 4 He, the spontaneously broken symmetry is a U(l) phase symmetry, 
in which the complex order field goes over into 

V>(x) - e^V(x). (2) 

In solids, the translational symmetry of the atomic position 

— > + u;(x), (3) 

is spontaneously broken. In ferromagnets, it is the rotational symmetry of the direction 
m of the elementary magnets 

m(x) — >■ i?m(x). (4) 

which is spontaneously broken. 

As the temperature is raised, the T D -behavior goes over into the typical Debye 
type of behavior. In crystals, it follows a universal function of T/0£>, where 0£> is the 
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following the Dulong-Petit law which depends only on the number of harmonic degrees 
of freedom. 

This picture is modified by anharmonic effects which become visible at increasing 
temperatures and which may be neglected in a qualitative discussion. In addition, 
there are usually massive excitations which, however, are strongly suppressed at low 
temperatures by a Boltzmann factor y = e~ E ^l T called the fugacity of the excitations. 
Here E a is the activation energy defined as the energy of an excitation of infinite wave 
length, which in relativistic systems is called the mass of an excitation. 

At a certain temperature T c , which in solids lies above the Debye temperature, in 
quantum solids and superfluid 4 He near or below the Debye temperature, the statistical 
behavior undergoes a dramatic change. The physical properties of the system change 
in a phase transition. This often happens at energies where the fugacities of massive 
excitations are very small so that up to the transition, only iVG-modes should be 
substantially excited. The dramatic change is caused by the sudden appearance of line- 
like excitations, whose large activation energies are overcome by a high configurational 
entropy. In solids, these excitations are defect lines, in superfluids, they are vortex 
lines. This entropy mechanism mechanism was noted in the fifties by Shockley and 
Feynman. A line-like excitation of length / with an excitation energy e per unit length 
is suppressed by a Boltzmann factor e~ d / T . If the line can bend easily on a length scale 
a, it has roughly (2D) 1 ^ possible configurations We can therefore estimate the partition 
function of an ensemble of such lines with various length as 



Above this temperature, the integral diverges and the ensemble undergoes a phase 
transition in which the lines proliferate and become infinitely long. We shall call this a 
condensation process of lines. 

It is an interesting task to construct a theory which explains the phase transitions 
in NG-sytems with the help of these line-like excitations without invoking explicitly the 
non-linear microscopic properties which are the origin of these excitations. The purpose 
of these lectures is to show that this is possible: There exists a minimal extension of the 
Debye theory of the harmonic NG-excitations which exhibit a phase transition caused 
by the condensation of idealized infinitely thin line-like excitations. The important the- 
oretical tool to describe these is a new type of fluctuating fields which I call generically 
defect gauge fields, in the superfluid vortex gauge field. In contrast to ordinary fluctuat- 
ing fields, vortex or defect gauge fields have the useful property that they are capable if 
explaining the above phase transitions with a quadratic field action. This makes it easy 
to transform their partition functions to various other forms by quadratic manipulation 
of functional integrals. By a particular variable change, the partition function can be 
converted into an interacting field theory of the |?/>| 4 -type which is called disorder field 
theory. 

The main difference with respect to the familiar order field theories of the Landau 
type is an opposite temperature behavior. In contrast to order fields, disorder fields 
acquire a nonzero expectation value at high temperatures. 




(5) 



The integral converges only 




(6) 
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to long-range interactions between the line elements. The interaction can be explained 
by a local field interaction if the description of the massless NG-modes is changed into 
a gauge field description. This is always possible. The new gauge fields are called 
generically stress gauge fields. The associated gauge invariance is dual to that of the 
defect gauge fields. 

In this way we arrive at a stress- and defect-gauge invariant theory of Goldstone 
modes and defect lines, or at an equivalent stress-gauge-invariant theory of disorder 
fields. 

Since NG-modes occur as a result of a spontaneous breakdown of some continuous 
symmetry, their complex order field description has the form 

V>(x) = e -WVo (7) 

where ip is a fixed order parameter and G a is a subset of generators of the associated 
symmetry group g. In the limit of infinite wave lengths, # a (x) becomes a constant 
and the energy is independent of 6 a — this being the reason for the masslessness the 
NG-modes. In any order field theory of the Landau type, the so-called London limit 
selects from all excitations precisely the NG-modes and approximates their field energy 
harmonically, while freezing all other excitations. The gradient energy of the order field 
gives rise to a gradient energy of the NG-modes in the London limit 

i NG oc J d D x\de-^ Ga \ 2 . (8) 

It is a nontrivial to translate this energy explicitly the # a (x)-fields. There always exist 
points in the field space a for which the order fields e - 8 M x ) G< * are identical. The field 
# a (x) can therefore jump between such values without any cost in gradient energy. If 
we want to describe the statistical mechanics of NG-bosons, these jumps have to be 
included into the partition. This will be done with the help of vortex or defect gauge 
fields as we are now going to explain with the help of several examples. 

The content of these lectures is largely contained in my textbook on this subject 
[1], except for discussions of some recent development. 



2. Superfluid 4 He 

The simplest system containing NG-modes is superfluid 4 He. They are describes by the 
phase factor of a complex order field 

V>(x) = e^Vo- (9) 
The line-like excitations are vortex lines. 

2.1. Gradient Energy 

In Landau's complex order field theory of superfluid 4 He, the invariance under U(l) 
transformation ^(x) — > e 8Qf ?/>(x) is spontaneously broken. As a consequence, the phase 
fluctuations of ^(x) are massless. On account of its exponential nature, the order field 
(9) has the following important topological property: For every point x it is possible 
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properly, let us introduce a ^-function on a volume V 

5(x; V) = / d 3 M (3) (x-x). (10) 



Then the above topological property of the order field amounts to an invariance under 
the replacement 

0(x) -»■ #(x) + 27r<5(x;V). (11) 

When expressing the gradient energy (8) in terms of 8$, this invariance must be main- 
tained. 

Moreover, the freedom of jumping by 2tt on a certain volume V opens up the 
possibility of #(x) jumping by 2tt across some higher-dimensional manifold, for instance 
a surface S in the three-dimensional superfluid. Due to the identity of the order field 
on both sides, the surface does not carry any energy. At the boundary L of the surface, 
however, the gradient energy de' 6 ^ does not behave smoothly. The boundary line 
carries an energy. It is a vortex line. To describe the energy of such a configuration we 
define a ^-function on a surface in D dimensions 

dxj(s, t) dxj(s, t) 



^ij(x; S) = J dsdt 
In three dimensions, we form 



ds dt 



st) 



6W(x-x(s,t)). (12) 



4(x; S) = -e ijk 8 jk (x; S). (13) 



1 
2 

It is now easy to see that the euclidean action associated with the correct gradient 
energy of the NG-modes in the presence of a vortex line around S must read 

Ang = ^tJ d3x{dl9 ~ 91)2 (14) 

where 

e^x) = 2irSi(x;S) (15) 
removes the singularities in <9;#(x) arising at the physically irrelevant jumps of the field 
0,-(x). 

Since the surface S is physically irrelevant, it can be deformed into any shape 
without changing the energy, as long as the boundary line stays the same. Indeed, 
under such a deformation S — > S' } the difference of the two surfaces encloses some 
volume V, and the ^-function (13) changes as follows: 

£(x; S) ^ £(x; S') = 5,-(x; S) + d,-5(x; V), (16) 

For the field #J(x), this implies 

0r(x)->0r(x) + 0,-A fi (x) (17) 

with the gauge function 

A fi (x) = 27r<5(x;V). (18) 
The energy (14) is invariant under this, if #(x) is simultaneously transformed as 

0(x) ^ 0(x) + A fi (x) (19) 

The field 0J is called the vortex gauge field of the superfluid and the combined trans- 
formation (17) and (19) is referred to as a vortex gauge transformation . 

The subscript 8 of the gauge function (18) emphasizes the fact that A§(x) is not 
an arbitrary function but a ^-functions on an arbitrary volume. 
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The physical content of the vortex gauge field #J(x) appears when forming its curl. By 
Stokes' theorem we find 

eijkdj8 k (x; S) = 6;(x; L) (20) 

where 

8i(x:L) = I «( 3 Vx-x) (21) 



is the ^-function on the boundary line L. The curl of 0j{x) is called the vortex density 

8 x v (x) = j v (x) = 2tt6(x; L) (22) 
As a consequence, the vortex density satisfies the obvious conservation law 

9-j v (x) = (23) 

This states that vortex lines are closed. 

The conservation law is a trivial consequence of j v being the curl of V . It is a 
Bianchi identity associated with the vortex gauge field structure. 

The expression (14) may not yet be the complete energy of a vortex configuration. 
It is possible to add a gradient energy in the vortex gauge field. It introduces an extra 
core energy for the vortex line. The extended action is 



*4ng = / d 3 x 



^(8tf-0 v ) 2 + |(8x0 v ) 2 



(24) 



The core energy term is a square of a ^-function. It must be regularized by using 
a slightly smeared-out ^-functions over some small length scale a. This length scale 
will be set equal to unity in the entire lecture. The regularizes last term is obviously 
proportional to the length of the lines. 

2.3. The Partition Function 

The partition function NG-modes and all possible vortex lines may be written as a 
functional integral 

Z NG = / D#]T e -^ G . (25) 

J {S} 

The measure J T>0 is defined, as usual by discretizing the space into a simple cubic 
lattice of spacing an integration at each lattice point over all 6 £ ( — oo, oo), and taking 
the continuum limit a — > 0. The sum over all surface configurations is defined on the 
lattice ba setting 

^(x;5') = 27rn t -(x) (26) 

on each lattice point, where n 8 (x) is an integer- valued field, and by summing over all 
integer numbers n 8 (x) 

E = £ • < 27 > 

{5} { ni ( x )} 
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latt 



E 

{ra 8 (x)} 



n / mm 

:e ~7 I <M|[9e(x)-2™(x)] 2 + f [Vxn(x)p}_ 



f28) 



which for e c = is the famous Villain approximation to the XY-model. The lattice 
spacing a has been set equal to unity and the symbol V denotes the lattice derivative. 

The vortex gauge fields requires a gauge fixing. On the lattice we can enforce the 
axial gauge 

n 3 (x) = 0. (29) 

The Lorentz gauge V • n = cannot be chosen for integer- valued fields. 

Similarly, the functional integral in the partition function (25) with the action (14) 
or (24) requires a gauge-fixing functional $[8 V ], for which we choose the ^-functional 



$[0 V ] = sm 



(30) 



Due to the sum over the vortex gauge fields V , the partition function (25) describes 
superfluid 4 He not only at zero temperature, where the NG-modes were identified, but 
at any not too large temperature. In particular, the phase transition is included. The 
vortex gauge field enriches the ensemble of fluctuating NG-modes in precisely the same 
way as the size of the order field ip does in a Landau description of the phase transition. 
In fact, it is easy to show that that near the transition, the partition function (25) can 
be transformes into a |?/>| 4 -theory of the Landau type. 

2.4. Interaction Energy between Vortices 

Let us calculate the interaction energy between vortex lines. Omitting the core energy, 
we expand the action (14) into 



*4ng = / d 3 x {{diOf + 268,6] + 6f 



(31) 



A quadratic completion leads to 



-4ng = ^ / d x 



di[9 + 



-d 2 



ch6J 



+ or - m- 



8 



2 d,6J 



(32) 



In the partition function, the 6 field can be integrated out producing a factor 
[Det(-8 2 )] ' . The remaining action can be written as 



A 



NG 



— / d 3 x(d x V )^ T (8 x V ^ 
2T J y ; -8 2V 

1 



2T J J -8 2J 



(33) 



This has the form of the magnetic Biot-Savart energy for current loops, except for an 
opposite sign (parallel currents attract each other, parallel vortices repel each other). 
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The gauge freedom in the partition function (25) is destroyed if the original U(l) sym- 
metry is broken. This adds a mass term to the NG-modes, and the action (14) becomes 



A 



NG 



2T 



d 3 x {[d,-0(x) - (x)f + m 2 #(x) 



(34) 



As a consequence of the gauge dependence, the previously irrelevant surfaces become 
energetic. To see this we write the action as 



A NG = ^J d 3 x { [(d t 9) 2 + m 2 9 2 ] - 9d t 9J + 9J 2 ] }. 



Integrating out the #-field in Z leads to an interaction 



•/4-int — _ m / d X 



This can be rewritten as 



Ant = / d 



1 

2T 



f8 x V ] 



iv2 



1 



1 



-8 + w? 



-8 + m 2 



(8 x V ) + m 2 9 



2 m 



1 



-8 + w? 



-9J 



(35) 



(36) 



(37) 



The first term is of the Biot-Savart type: antiparallel line elements attract each other 
with a finite-range Yukawa potential e~ mr /Airr. The second term describes a Yukawa 
interaction between the normal vector of the surface elements. It produces a physical 
surface of thickness 1/ra. 

The same gauge-breaking mechanism can be used to construct a simple model of 
quark confinement which will be described in Section 4. 



2.6. Gauge Field of Superflow 

There exists an alternative representation of the partition function in which the vortex 
lines appear directly via their physical vortex density; the jumping surfaces S being 
completely eliminated. This is possible by representing the NG-modes in terms of a new 
gauge field. It is canonically conjugate to the angular field 9 and called generically the 
stress gauge field, which here is the gauge field of superflow . Recall that the momentum 
variable in an ordinary path integral 



Vxe-f dt 



(38) 



is introduced by a quadratic completion, rewriting (38) as 



DxDpe 



fdt 



tpx- 



By analogy, we rewrite the partition function (25) as 

Z = j V9 j Vbi e ~ Al 



(39) 



(40) 



{S} 
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^ng = J d 3 x |-6 2 (x) - ibi [0,-0(x) - ^ v (x)] j • (41) 

Note that if we go over to a Minkowski space formulation in which 
plays the role of the time, the integral 

creates on a discretized time product of ^-functions 

(On+l \0n)(0n\0n-l)(0n-l l^n-2) (43) 

with 

(e n \e n _ 1 ) = s n (e n -e n _ 1 ) (44) 

which can be interpreted as Dirac scalar products in the Hilbert space of the system. 
On this Hilbert space, there exists an operator 6;(x) whose zeroth component is given 
by 

b = -idg (45) 
and satisfies the equal-time canonical communication rule 

6o(xi, x ), 0(x'_l, so)] = ^ (2) (x ± - x'J. (46) 

The charge associated with &o(x), 

Q(x ) = j d 2 xb (x. L ,x ), (47) 

generates a constant shift in 6: 

e -«-Q(x )^( x _ Lj Xo ) e i°>Q(.*>) = #( x± , x ) + a. (48) 

Thus it multiplies the original field e ld ^ by a phase factor e ia . The charge Q( x o) 
is the generator of the U(l) symmetry transformation whose spontaneous symmetry 
breakdown is responsible for the NG-nature of the phase fluctuations #(x). Since the 
original theory is invariant under the transformations ip — > e ta if) } the energy (41) does 
not depend on 6 itself, only on d$. Thus, in the partition function, the field 6 can be 
integrated out using the formula 

' VejW = 8[f(x)]. (49) 
This produces the conservation law 

<%(x) = 0, (50) 

implying that Q( x o) is a time independent charge and e ta ® a symmetry transformation. 

If the energy does depend on 6, the charge Q( x o) becomes x -dependent. However, 
it still generates the above U(l) transformation. 



1 Similar canonical representations for the defect ensembles in a variety of physical systems in Ref. 

[2]- 
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representation (41) contains more complicated functions of 6;(x). 

The conjugate variable to the phase angle of a complex field is the particle current, 
which is here the particle current in the superfluid condensate. Thus 



j,-(x) = 6 t -(x) 



(51) 



can be interpreted as the supercurrent . 

The dynamical conservation law <9;&;(x) = can be ensured automatically as a 
Bianchi identity, if we represent b curl of a gauge field of superflow 



6,-(x) 

Then the action (41) goes over into 



ijkdjcikfe). 



NG 



d x 



- v 8xa) 2 + za.(8x0 v ) + ^(8x0 v ) 2 



The partition function becomes 



Z= I Da]Te-^ 

{5} 



(52) 



(53) 



(54) 



The action is doubly gauge invariant, under the gauge transformations of superflow 

a,-(x) -»■ a,-(x) + d,-A s (x). (55) 

and under the vortex gauge transformation (17) and (18): 

0y(x)->0r(x) + 0.-A fi (x) (56) 

with the gauge functions 

A fi (x) = 2ttS(x;V). (57) 

With respect to the first gauge structure, the current conservation law <9;&;(x) is a 
Bianchi identity. 

The integrations in the partition function (54) require a gauge fixing factor for 
both gauge fields. For the vortex gauge field, we use the axial factor (30), for the gauge 
field of superflow, any fixing will do, for instance the transverse one: 



$[a] = 6[d ■ a]. 



The action A% G can also be written short as 



•A. — I d x 



|(8xa) 2 + za.r + ^r 2 



(58) 



(59) 



In this expression, the jumping surfaces have disappeared and the action depends only 
on the vortex lines. The partition function may be rewritten as 



Z = I Da]T e -^ 

J {L} 



(60) 
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$[n = *[8.j v ] (6i) 

which ensures the closedness of the vortex lines. 

For a fixed set of lines L, the partition function (60) has a similar form as the 
magnetic partition function of a set of current lines. Around a vortex line, the field 
b(x) = 8 X a(x) looks precisely like a magnetic field B(x) = 8 X A(x) around a current 
line. Thus, integrating out the vector potential a yields once more the same Biot-Savart 
interaction energies as obtained before in another way in Eq. (33). 

Note that if the energy contains an explicit dependence like the mass term of the 
last section, there exists no a-field representation of the NG-modes. 



2.7. Disorder Field Theory 



The sum over all vortex line configurations in the partition function (60) can be done 
separately, defining an a-dependent vortex partition function 



ZM = ]T««exp 



{S} 



d 3 x 



2T 



8 x V V + za- (8 x V ^ 



J^<5[8 • j v ] exp 

{L} 



d 3 Q 



c c . v 2 . ■ -v 

T^J + «a • j 



(62) 



This can be viewed as a particular representation of an entirely different set of NG- 
modes, the dual NG-modes. To see this we make use of an important relation: 



{L} {L} 

This can easily be proved by going on a lattice and using the Poisson formula 

n m 

It is therefore possible to rewrite the partition function (62) as follows 

d 3 x 



(63) 



(64) 



Z v [a] = J2 Vbl / ^ ex P 
{L} J 

Integrating out the 6 8 (x)-field gives 
Z v [a] = J2 I ^ ex P 



^jV + iaA + —jj ( m - ej 



(65) 



1 T 



(66) 



This is the NG-theory of another [/(l)-invariant |^) 4 |-theory, the above-mentioned dis- 
order field theory. It differs from the initial |?/>| 4 theory by the presence of an external 
vortex gauge field a in the gradient term coupled via a covariant derivative. The disorder 
field partition function is 



Z v [a] 



DipDip* exp 



d 3 x 



l/fl o \ |2 , 2 I |2 . 9 I |4 

\{Oi - Zira,) tf\ + m \<p\ + - \<p\ 



(67) 
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frozen out. The functional integral over 6 removes the longitudinal part of a 4 -, and we 
remain with 



Z v [a] = exp J d 3 ia^j (68) 

where ^ 

m 2 a = -. (69) 

The exponent in (68) is a mass term for the transverse projection of the gauge field of 
superflow 

aj = - didi/Q 2 ) a r (70) 
In the complex disorder field theory (67), the mass term m 2 is negative and the disorder 



field acquires a nonzero expectation value ip = y'—m 2 /g. This produces again a mass 
term (68). 

The mass term (68) shows that in the high-temperature phase, the singular nature 
of the vortex gauge field V is lost and V can be treated as an ordinary fluctuating 
transverse field. The sum turns into a functional integral over an ordinary gauge 

field J X)0 V . Indeed, when doing this integral in the partition function out (62), we find 
that the vortex density has the simple correlation function 

< j7(x)j/(x') >= -(% - <W8 2 )«5( 3 )(x - x'), (71) 

which produces directly the mass term (68). 

This is to be contrasted with the contributions to the a-field action coming from 
Z v [&] of Eq. (62) at low temperatures. There the vortex lines form small loops, and 
these can be shown to give 2 

Z v [a] ~ exp |-i(8 x a)/(-z8)(8 x a) j . (72) 

Such a contribution changes only the dispersion of the gauge fields of superflow. In- 
finitely long vortex lines are necessary to change the mass. These are present at high 
temperatures where the correlation function of the vortex densities is (71). 

With the help of the disorder partition function Z v [a], the partition function (25) 
of the NG-modes associated with the action (24) can be replaced by another completely 
equivalent one: 

z NG = [v^Y! j vee ~^ G ( 73 ) 

{L} 

with the action 

A NG = T j d 3 x |i(8 x a) 2 + || (d6 - V - 2vra) 2 j . (74) 

This action is invariant under the following two types of gauge transformation. First, 
there is invariance under the gauge transformations of superflow (55), if it is accompa- 
nied by a compensating transformation of the angular field 6: 

a,-(x) -> a,-(x) + d,-A(x), 0(x) -> 0(x) + 2ttA(x). (75) 



2 See Ref. [1], Vol. I, Part 2, Section 9.7 and Ref. [3]. 
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field, the analogs of Eqs. (17)— (19): 




^ v (x) + d,-A fi (x), 0(x) ^ 0(x) + A 5 (x), 



(76) 



with gauge functions 



A 5 (x) = 27r<5(x; V). 



(77) 



At high temperatures, the vortex lines in V are frozen out and the action (74) shows 
again the mass term (68). 

The mass term implies that at high temperatures, the gauge field of superflow 
possesses a finite range. At some critical temperature superfluidity has been destroyed. 
This is an analog of the Meissner effect within the disorder theory [4]. Note that in 
the absence of the gauge field of superflow a, the field 6 would be of long range, i.e., 
massless. The gauge field of superflow absorbs the massless mode and the system 
has only short-range excitations. More precisely, it can be shown that any correlation 
function involving local gauge-invariant observable quantities must be of short range in 
high-temperature phase. 

Take, for instance, the local gauge-invariant current operator of the disorder field 



This has no zero-mass pole. 

If we remove the locality restriction, it is easy to find a filed which has long-range 
correlations in the superconducting phase. If there were no magnetic interactions, the 
exponential filed e 10 would provides a order parameter of the system (which would then 



invariant. Its expectation value vanishes identically on account of Elitzur's theorem [5]. 
There exists, however, a simple gauge-invariant modification of this used first by Dirac 



j s = 89 - a 




In the gauge 6 = 0, we easily find 




(79) 





[6, 7]. Let 





be a charge density localized at the point x. Then the disorder field 





is gauge invariant and the expectation value 



6 = (elf x) ) 




is nonzero. 



Note that apart from being invariant under magnetic gauge transformations 
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At low temperatures where the vortex lines can be ignored, it can easily be shown 
to be nonzero. We calculate the Green function associated with two fields (82): 

^( Xl ,x 2 )^( e 2 (xi) e / (X2) ). (85) 

In the transverse gauge, we can easily integrate out the ^-fluctuations in (73) and find 
the same correlation function as in (111) 

G f o(xi,x 2 ) — > const. (86) 

|X1-X 2 |-KX> 

Hence 

6 = (e^) = const. (87) 

This is nonzero, just as the order parameter of the superfluid in the low temperature 
phase. 

2.8. Disorder Theory of Superconductor 

A superconductor is described by the Ginzburg-Landau action 



A G l = 4 J d 3 x { |(8 - zA) ^ + m 2 M 2 + 9 - M 4 + 1(8 x A) 2 } . 



f88) 



We have set the charge of the Cooper pairs equal to unity, for simplicity. In the 
superconducting phase, m 2 is negative. Going to the London limit, the field <p is 
approximated by 

<p = e i§ W<p (89) 



with a constant real Lp = y 7 —m 2 / 'g, and the action becomes 

1 



^(8^-e V -A) 2 + i(8xA) 



(90) 



where we have set (p^ = m 2 i /2e 2 . The partition function reads 

z sc = j ve j vA^2 e ~ Ahh - ( 91 ) 

The action (90) has the same form as the action in the disorder representation (74) 
of superfluid 4 He. The role of the gauge field of superflow is now played by the vector 
potential A of magnetism. The action has the following two types of gauge invariances: 
the magnetic ones 

A 8 (x) ^ A 8 (x) + d,-A(x), 0(x) ^ 0(x) + A(x), (92) 
and the vortex gauge transformations 

#7(x) -^r(x) + d,-A fi (x), 0(x)^0(x) + A fi (x), (93) 
with gauge functions 

A fi (x) = 2ttS(x;V). (94) 
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function (91) describes the statisctical behavior of the superconductor not only at zero 
temperature, where the action (90) was derived, but at all not too large temperatures. 
The fluctuating vortex gauge field V ensures the validity through the phase transition. 

We shall now derive the disorder representation of this partition function in which 
the vortex lines of the superconductor play a central role in describing the phase tran- 
sition. 

At low temperatures, the vortices are frozen, and the ^-fluctuations in the partition 
function can be integrated out. This reduces the action to 



Al ~ ^ J d 3 xA 2 T , (95) 

i.e., to a simple transverse mass term for the vector potential A. This is the famous 
Meissner effect in a superconductor, which limits the range of a magnetic field to a finite 
penetration depth X = 1/ra^. The effect is completely analogous to the one observed 
previously in the disorder description of the superfluid where the superfluid acquired a 
finite range in the normal phase. 

To derive the disorder theory of the partition function (91), we supplement the 
action by a core energy of the vortex lines 



Acre = ^ j d 3 x(d X )\ (96) 

As in the partition function (40), an auxiliary hi field can be introduced to bring the 
action (90) to the form 



All = / d 3 ; 



-^b 2 + ib (d 9 - 6 V - A) + -i (8 x A) 2 + -^(8 x V 

2m 2 , V / 2T 2T 



(97) 



By integrating out the ^-fields, we obtain the conservation law 

8 • b = (98) 
which we fulfill by introducing the gauge field of superflow in the superconductor 

b = 8 x a. (99) 



In this way, we find the action 

f _ .„ _ 1 



r-3 

All = I x 

where 



, ,8 x a) 2 — za • (8 x A) + -4^(8 x A) 2 - za • j v + 

2m 2 / ; V ' 2T ' 2T 



(100) 



j v = 8 x . (101) 

At low temperatures T where vortex lines are absent, the a-field can be integrated out 
and yields once more the transverse mass term (95). 

At high temperatures T, the vortex lines are prolific and the vector vortex gauge 
field V can be integrated out like an ordinary field. With the help of the correlation 
function (71), we obtain the transverse mass term 

T f m 2 

d 3 x^a 2 T (102) 



2m\ 
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(103) 



When integrating out the massive a-field, no mass term is produced for the vector 
potential A. Instead, the vector potential acquires the effective action 

A A = -i / d 3 x 8 x A (\ + ^ ^ 8 x A. (104) 

2T J V -B 2 + ml) y ' 

Expanding the denominator in powers of — 8 2 we see that only gradient energies appear. 
Thus, apart from a modification of the dispersion, the vector potential is of long range 
giving Coulomb-like forces at large distances. 

For the low-temperature mass value m 2 = 0, the second term reduces again a 
transverse mass term (95) associated with the Meissner effect. 

We can represent the fluctuating vortices in the superconductor by a disorder field 
theory in the same way as we did for the vortices in the superfluid, i.e., we repeat the 
transformations from (62) to (66). The angular field variable of disorder will be called 
6, the vortex lines in the disorder theory V . The disorder action reads 



T (8 x a) 2 — z'a • (8 x A) + i(8 x A) 2 + f m 2 (8fl - V - 2vra) 2 



2m\ 2T 8ir 2 m 2 A 

Near the phase transition, this is equivalent to a disorder field action 



(105) 



^ll~ /^<|^(8xa) 2 -za.(8xA) + -l(8 
' 1 2m A 2T 



+ |(8 - 2vra) Vf + m 2 |Vf + |H 4 } • (106) 

The loop diagrams of the disorder field ip are the pictures of the vortex lines. 

Note that on the average, za is equal to the magnetic field 8 X A. 

The disorder fiels theory for the superconductor has recently led to a determination 
of the critical properties of the superconductive phase transition in a 4 — e expansion 
[8]. Within the Ginzburg-Landau theory, this is impossible since the renormalization 
group gives no infrared stable fixed point at all. 

2.9. Order versus Disorder 

Superfluid 4 He 

The order parameter O of the original theory (14) of (24) is defined as the expectation 
value of the phase of the complex order field 0(x) = e ld ^\ 

O = (0(x)) = (e^ x )). (107) 

Its behavior can best be studied by considering the correlation function of two order 
fields O(x): 

Go(x,x') = (c'^c^M). (108) 

In the low temperature phase, where vortices are rare the #(x)-field fluctuates almost 
harmonically. By Wick's theorem, we can then approximate 

G (x 1 ,x 2 )^ e -^K e ( Xl )- e ( X2 )] 2 >. (109) 
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goes to zero for r — s- oo. After a renormalization of the fields, we see that 

Go(xi,x 2 ) — > const. 

|xi-x 2 |-kx> 

By the cluster properties of correlation functions, 

(0( Xl )0(x 2 )) — ► (0(x 1 ))(0(x 2 )), 

|Xl-X 2 |-KX> 

this implies that the order parameter is a constant: 

O = (O(x)) = const. 



liO) 



111] 



112) 



113) 



Let us calculate the order parameter in the high temperature phase. To find the 
correlation function Go(xi,x 2 ), we insert into the partition function a term 



e -t 3 (x)0(x) 

with two point sources of opposite charges at the points Xi and x 2 : 

g (x) = «5( 3 )(x-x 1 )-«5( 3 )(x-x 2 ). 
When going to the canonical representation (41), we obtain an action 



A= / d 



| 6 2_ i6 . (a . 6) _ 6 ,v ) + i9(x)6)(x) 



Integrating out the 6 field in the partition function gives the constraint 

<9A(x) = -q(x). 



114) 



115) 



116) 



117) 



In the magnetic-field analogy, the point charges in <?(x) correspond to a magnetic mo- 
nopole and an anti-monopole. 

The constraint is solved by introducing a singular field 

&T(x) = -Mx;Z) (118) 

where the line L runs from Xi to x 2 along an arbitrary path. This field satisfies 

d^(x) = - 9 (x). (119) 

We can therefore write 



b(x) = 8 x a(x) - b m (x) 



and the action (116) can be replaced by 



A= d 



- (8 x a - b 11 

2 v 



\2 ■ -v i e c . v2 



(120) 



121) 
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fixed ends results in a transformation 



&T(x) ^ &TM + e ijk dMx; s), (122) 

where S is the surface over which the line L has swept. The action (121) remains 
invariant if the gauge field of superflow is simultaneously transformes as 

a 8 (x)^a 8 (x) + ^(x;5'). (123) 

The pure phase order field e ld ^ corresponds to a gauge field 

6-(x) = ^(x;X 1 ) (124) 

with a line L\ starting at x and running to infinity, i.e., it is represented in the path 
integral (73) by the integrand 

e «*(x)^ e -T / d^{|b m (x) 2 -b m (x).[3xa(x)]}_ / 125 x 



When integrating out the gauge field of superflow a at low temperatures where 
the vortex lines are frozen out, the action the partition function contains b m in the 
form of a factor 



T 



2 jV^ bm2 _ (9xbm) _i_ (9xbm) j ^ ^_Zf d 3 xd . h ™_L_ d . h ™ 



From this we obtain the correlation function 

G (x 1 ,x 2 ) = e - f/d3 ^ 9 . (127) 

Inserting (119), this is precisely the correlation function (109) with (110), which goes 
to a constant for |xi — x' 2 | — > oo as in (111) implying a nonzero order parameter. 

In the disorder version Zng of the partition function Zng, it is easy to calculate 
the correlation function also in the high-temperature phase. There, the a-field contains 
the transverse mass term (68) and the situation is very similar to that in Section 2.5. 
The line connecting the points Xi and x 2 becomes physical and its energy grows linearly 
with the length of the line. Thus, the correlation function vanishes 

G (Xl,X 2 ) ~ e -c°n S t.| Xl -x 2 | Q. ( 12 g) 

|X1-X 2 |-KX> 

By the cluster properties (112) of correlation functions, this shows that the order pa- 
rameter vanishes. 

Note that precisely the same exponential falloff is found within Landau's complex 
order field theory: 

< V(xr)V(x 2 ) >oc / jf* = -l ^'"""' . (129) 

J (Ztt) 6 k 2 + m 2 47r |xi — x 2 | 

The mechanism by which the finite range arises is, however, quite different. The size 
fluctuations of the order field play an essential role. In the partition function (25), their 



3 See Ref. [1], Vol. I, Part 2, Chapter 10. 
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end of Section 2.3. 



b) Superconductor 



The order parameter of the superfluid supplies directly a disorder parameter of the 
superconductor [9]. Inserting the monopole gauge filed (124) into the actions (100), 
(105) and (105), we obtain the disorder field analogous to (125): 

£)(*) = e -?/ d3 4l bm ( x ) 2 - bm ( x )-[ 9xA ( x )]}. (130) 

The action for the correlation function of a monopole antimonopole pair is 



A 



d 3 x 



5-(8 x a) 2 — za • (8 x A) H (8 x A - b m ) 2 - za • j v + 

2m\ K ' K ' 2T ' 2T ° 



31) 

with b m of Eq. (118). At low temperature T of the superconductor, where the vortex 
lines of the superconductor are frozen out and the integral over the a-filed yields a mass 
term for the vector potential A, this becomes 



A 



d A x { -L 
2T 



(8 x A - b m r + -rp Ay 



za • j H -r 



(132) 



Just as in (128), the line connecting the endpoints acquires an energy and the correlation 
function has the property 



G i5 (x 1 ,x 2 ) ~ e- const - |xi - X21 — > 0. 

|X1-X 2 |-KX> 

implying a vanishing disorder parameter: 

V = (D(x)) = 0. 



(133) 



(134) 



In the normal phase of the superconductor, the monopole gauge field enters the 
action (104) which becomes 



A! \ 



1 

2T 



d 3 x (d x A - b r 



1 + 



m 



-8 + m\ 



f8 x A — b 11 



(135) 



As in (127), we obtain for a large distance between the monopole and the antimonopole 
a Coulomb interaction implying that the disorder parameter is a constant: 



V = (-D(x)) = const. 



(136) 



2.10. Order of Superconducting Phase Transition — Tricritical Point 

In the superconductive phase, there are only a few vortex lines and the field ip fluctuates 
around zero. The mass term is m 2 > 0. The vector potential A may be integrated out 
and leads to an action for the a-field 



Aa 



dh 



—^(8 x a) 2 + -al + ai-IVH 2 
2m 2 / J 2 T T 2 m 



(137) 
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between a and ip simple. This yields also a coupling term a^|?/>| 2 involving the lon- 
gitudinal components of the field a. When integrating out the longitudinal part, the 
functional determinant is canceled by the gauge fixing factor of the gauge ?/>=real. 

The a-field can now be integrated out; for a constant ip they yield an effective 
action 



This term lowers the (<?/4) \t/j\ 4 interaction term in (106) . An increase in itia corresponds 
to a decrease of the penetration depth in the superconductor, i.e. to materials moving 
towards the type- 1 regime. At some place, the |?/>| 4 -term will vanish and the disorder 
field theory requires a |?/>| 4 -term to stabilize the fluctuations of the vortex lines. In such 
materials, the superconductive phase transition turns from second to first order. The 
existence and position of such a tricriticsl point was first established in [10]. 

2.11. Vortex Lattices 

The model action (14) represents the gradient energy in superfluid 4 He correctly only 
in the long- wavelength limit. The neutron scattering data yield the energy spectrum 
to = S(k) shown in Fig. 3. To account for this, the action should be taken as follows: 





When expanded in powers of \ip 



2 , this produces a term 




(139) 
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Figure 3. The structure factor of superfluid 4 He measured by neutron scat 
tering showing the excitation energy of the NG-bosons. 
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Figure 4. The phase diagram of a two-dimensional layer of superfluid 4 He. 
At higher fugacity, the vortices first condense to a lattice, which undergoes 
the Kosterlitz-Thouless vortex unbinding transition only after a melting tran- 
sition. 
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Figure 5. The experimental phase diagram of a two-dimensional layer of 
superfluid 4 He diluted with 3 He which decreases the fugacity and separates 
the vortex melting transition from the Kosterlitz-Thouless transition. 
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A NG = —j d 3 x(d 9 - e v )^^(8 9 - V ). (140) 

The roton peak near 20A" 1 gives rise to a repulsion between opposite vortex line ele- 
ments at the corresponding distance. If a layer of superfluid 4 He is diluted with 3 He, 
the core energy of the vortices decreases, the fugacity y and the average vortex number 
increases. For a sufficiently high average spacing, a vortex lattice forms. In this regime, 
the superfluid has three transition when going from zero temperature zero to the nor- 
mal phase, first condensation process of a vortex lattice, then a melting transition of 
this lattice into a fluid of bound vortex-antivortex pairs, and third a pair-unbinding 
transitions of the Kosterlitz-Thouless type [11, 12]. The latter two transitions have 
apparently been seen experimentally [13]. 



3. Crystals 

In crystals, the elastic energy is usually expressed in terms of a material displacement 
field u;(x) as 

E = j d 3 xul(x) (141) 
where ^ 



u 



t3y x) = -(d iUj (x) + d jUj (x)) (142) 

is the strain tensor. For simplicity, we have set the shear elastic constant // equal to 
unity and neglected the second elastic term A(<9;u 8 ) 2 . The elastic energy goes to zero 
for infinity wave length since in this limit it 4 -(x) reduces to a pure translation and the 
energy of the system is translationally invariant. The crystallization process causes a 
spontaneous breakdown of the translational symmetry of the system elastic distortions 
are iVG-modes of this symmetry breakdown. 

A crystallize material always contains defects. In their presence, the elastic energy 

" . «-/" -<> 

where it?- is the so-called plastic strain tensor describing the defects. Topologically, the 
crystal defects arise as a consequence of the multivaluedness of the displacement field 
on a lattice. Due to the identity of the lattice constituents, the displacement is defined 
only modulo lattice vectors. For mathematical simplicity, let us assume the lattice to 
be simple cubic of spacing 2tt. Then the energy of it 4 -(x) and it 4 -(x) + 27rA^-(x) must be 
indistinguishable for integer valued fields A^-(x) which correspond to a permutation of 
the lattice sites. The order field of a crystal has therefore the three components 

^■(x) = e MX Vo 8 . (144) 

This is very similar to the order field of a superfluid. There are now three types of 
surfaces where the distortion tensor fin = d{U\ jumps by 2tt } one for every lattice 
direction h\ = 1,2,3. They are describes by the plastic distortion 

Pf l (x) = 8 i (x;S)b l . (145) 

The surfaces S are irrelevant and correspond to a layer of atoms missing from a crystal. 
Only the boundary line is physical. They are called dislocation lines of Burger vector 
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dislocation lines are extracted from by forming the dislocation density 

«4 X ) = eimndmPnli*) = $i{L)h- (146) 

Being a curl, it satisfies the conservation law 

diau = (147) 

trivially as a Bianchi identity. The conservation law states that dislocation lines are 
closed. 

A closed Volterra surface is associated with a translation 

u 8 (x) -> u 8 (x) + 2ttM(x; V) (148) 

where V is the volume enclosed by the closed Volterra surface S. Under such an 
operation, the order parameters e lu '^ are invariant. 

In a real crystal, dislocation lines are able to stack up on top of each other forming 
line-like defects called disclination lines . Their geometric properties are most easily 
understood by noting that instead of the trivial translation (148) across a closed Volterra 
surface we can form a combined translation rotation plus 

u 8 (x) ^ [k + i?^-(x)] 5(x; V) (149) 

where is a discrete rotation under which the crystal is invariant (in a simple cubic 
lattice, Rij consists of products of rotations by 7r/2). The nonabelian nature of the 
trivial gauge transformations (149) makes a full discussion of the plastic deformations 
quite complicated. A useful approximation consists in assuming to be a small 
rotation and to reexpress (149) in the infinitesimal form 

u;(x) ->■ u;(x) + (k + e lqn tt q x r )8(x; V). (150) 

Assuming the presence of a closed Volterra surface, we extract the jumps in the distor- 
tion field 

fiki = id k ui(*x) = 4(x; S)(k + ei qr tt q x r ) - 6(x; V)ei qk Vl q . (151) 

The first term is defined also for an open surface in which case it becomes the plastic 
distortion of a general defect line 

ft v kl = 4(x; S)(k + ei pr tt q x r ). (152) 

By forming the curl of the distortion field (151), the term <5(x; V) drops out and we find 

OLii = e l3k d 3 d k ui(y.) = <5 8 (x; L)(bi + ei qr Sl q x r ). (153) 

This is called the dislocation density of the defect line. The rotation field 

Wj = ^e jk id u ui (154) 

has the following derivative 

d n wj = ^tjkidkd n ui = ^tjkidkPli + <5 n (x; S)Vlj. (155) 
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4>i 3 = 6 n (s)tt r (156) 
The entire tensor (155) is called the plastic bend-twist 

K lj = d -^j = ^tikidnPki + €y (157) 

The curl ol the plastic rotation is the disclination density 

On = ^jkdj(/) p kl = 6;(x; L)Vli. (158) 

Being a curl, this satisfies the conservation law 

diOu = (159) 

implying that disclination lines are closed. This is again a Bianchi identity of the defect 
system. The conservation law for dislocations is modified by the presence of disclination 
as follows: 

3jC^ik £imn@mn' (160) 

Dislocation lines can now end in disclination lines. 

An important geometric quantity characterizing dislocation and disclination lines 
is the incompatibility or defect density 

r]ij(x.) = e lk ie jmn d k d m ui n (x). (161) 

It can be decomposed into disclination and dislocation density as follows: 

?/ij(x) = %(x) + ^<9 m [e min a jn (x.) + (ij) - e ijn a mn (x.)] . (162) 

This tensor is symmetric and conserved 

<9 8 %(x) = 0, (163) 

again a Bianchi identity of the defect system. 

The tensors ctij,6ij and rjij are linearized versions of important tensors in a 
Riemann-Cartan space , a noneuclidean space with curvature and torsion. Such a 
space can be generated from a flat space by a plastic distortion, which is mathemati- 
cally represented by a nonholonomic mapping 



Zi + u,-(x). (164) 



Such a mapping is nonintegrable. The displacement fields and their first derivatives fail 
to satisfy the Schwarz integrability criterion: 

(didj - djdi) u(x) ^ 
(didj-djdjdkmix) + 0. (165) 

The affine connection of the geometry in the plastically distorted space is 

Fiji = didjui. (166) 
24 



torsion tensor being 

Skij = (Tf/fc — Tjik)/2. (167) 
The dislocation density ctij is equal to 

ctij = tikiSkij, (168) 

The noncommutativity ol the derivatives in Iront ol <9fcU/(x) imply a nonzero cur- 
vature. The disclination density 6ij is the Einstein tensor 

9 %3 = G %3 = 3% -^g ]t R (169) 

ol this geometry. The tensor rji 1} finally, is the Belinlante symmetric energy momentum 
tensor which coincides with the Einstein tensor obtained Irom the Christoffel symbols 
ol the metric 

9ij = % + ®i u i + fy u i (I 70 ) 

rather than Irom the alfine connection To verily this statement we just note that 

the spin current density ol the gravitational field is (setting the gravitational constant 
equal to unity) 

= —2 (Sijk + SikSj — 8 3 kSi) . (171) 

For more details on the geometric aspects see Vol. II ol [1], Part IV, where the lull 
one-to-one correspondence between delect systems an Riemann-Cartan geometry is de- 
veloped as well as a gravitational theory based on this analogy. 

It is possible to write down an elastic energy which disentangles dislocations and 
disclinations by including higher gradients ol the displacements field. This energy reads 



E = d 



J3, 



(u l3 - ulf + e (d^ - kI) 2 \ . (172) 

The parameter i is the length scale over which the crystal is rotationally stiff. 

The partition lunction contains integrals over it 4 - and sums over the jumping sur- 
laces ol dislocations and disclinations. By integrating out the u 4 -fields, one obtains a 
Biot-Savart type ol interaction energies between the delect lines in which dislocation 
line elements interacts with each other via a Coulomb potential, and disclination line 
elements via an r-potential. 

It is again possible to eliminate the jumping surfaces Irom the partition lunction 
by introducing conjugate variables and associated stress gauge fields. For this we rewrite 
the elastic action ol delect lines as 



E = d 3 x 



Li , ^1 

L73) 

and form the partition lunction by integrating over cr 8J , r 8J , Ui^Uj and summing over all 
jumping surfaces S in the plastic fields. A lunctional integral over the antisymmetric 
part ol (Tj-j has been introduce the obtain an independent integral over Ui (if we were 
to integrate out the antisymmetric part of cr 8J , we would enforce the relation Ui = 
^Qjk{djUk + flfj)- By integrating out ujj and Ui, we find the conservation laws 

di(Tij = -tjkim, diTij = 0. (174) 
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and (159). 

They can be enforced as Bianchi identities by introducing the stress gauge fields 
Aij and hij and writing 



tikidkAij 
tikidkhij - 



j Au 



Aji. 



(175) 



This allows to reexpress the energy as 



E 



d 3 x 



(176) 



The stress gauge fields couple locally to the defect densities which are singular on the 
boundary lines of the Volterra surfaces. In the limit of a vanishing length scale £ } is 
forced to be identically zero and (175) allows us to express Aij in terms of hij. Then 
the energy becomes 



E 



1 

- (cr 

4 



((Tij + CTjif + hijljij 



(177) 



where the defect density rjij contains dislocation and disclination lines. 

Depending on the length parameter £ of rotational stiffness, the defect system was 
predicted to have either a single first-order transition (for small £), of two successive 
continuous melting transitions. In the first transitions, dislocation lines proliferate and 
destroy the translational order, in the second transition, disclination lines proliferate 
and destroy the rotational order [4]. 

In two dimensions, the existence of such a scenario (where the continuous transi- 
tions for small £ would be of the Kosterlitz-Thouless type) was conjectured a long time 
ago [14, 15], but the defect models constructed to illustrate this displayed only a single 
first-order transition [16]. The computer simulations of the model containing the length 
parameter £ [17], on the other hand, confirmed the theoretical expectation citekl2tr. 



4. Abelian Quark Confinement 

In Section II we have observed the appearance of confining forces caused by the break- 
down of vortex gauge invariance: An irrelevant geometric line becomes energetic and 
creates a linearly rising potential between the endpoints. This mechanism can be used 
in four spacetime dimensions to construct an abelian theory of quark confinement. In a 
superconductor, whose ground state contains a condensate of electric charges, magnetic 
monopoles are confined. To obtain confining force between electric charges, we have to 
construct a dual theory containing poliferating monopole world lines. 
The euclidean field action of such world lines is given by 

A=\j d*x{ [i^(x) - F-(x)] 2 + ij^x^x)} (178) 

where = d ll A v — d^A^ is the usual electromagnetic field tensor, 

j ll (x) = eS i (x;L) (179) 
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Figure 6. Separation of first-order melting transition into two successive 
Kosterlitz-Thouless transitions when increasing the length scale £ of rotational 
stiffness of the defect model. 




Figure 7. The phase diagram in the T-£-plane. 
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is the gauge field of mono-poles . The ^-function on the surface S 



<5a«(x; S) 
has the property that 



dodr 



3x m dx v , , 

"a » ^ 

do ot 



^ 4 )(X-X((7,T)) 



<9 k <5 Ak (x; £) = <5 A (x; L) 
where L is the boundary line of S. Hence 

1 



e^xJlF^ = ^(x) 



180) 

;i8i) 
;i82) 

;i83) 
;i84) 



where 

j'^x) = gS^(x;L) 

is the magnetic current of the monopole world line along L. The field energy invariant 
under monopole gauge transformations S — > S 1 under which 



where 



<5 Ak (x; S) ->■ <5 Ak (x; S) + e AK(7T 9 (7 (5 T (x; V) 



Qx Qx d X 

<5 T (x; V) = e^ K1 I do dr A— ^ — ^ — 6 {n] '(x - x(<r, r, A)) 



185) 



186) 



da dr d\ 

is the ^-function on the volume over which S has swept. The vortex gauge field trans 
forms as follows 

F™ - F£ + d,A™ - d„A 

where 



kf = ^ q S(x;V). 
The gauge functions A™ can be absorbed into the vector potential 

— > + k™ . 



187) 
188) 



189) 



The additional interaction term J J 4 xj |U (x)A™(x) changes by eg times the integral 
J d 4 xS ll (x; Z)(5 |U (x; V). The integral is an integer number n counting how many times 
the line L pierces the volume V. Thus, the functional integral of the system 



Z 



-A 



(190) 



picks up a phase factor e ie9 ' a . This is an irrelevant factor if 

eg = 27m, 
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which is Dirac's charge quantization condition. 4 By introducing an auxiliary field j ^ 
the theory can be transformed to the dual form 



A 



d 4 Q 



(192) 



4 Dirac uses a field energy with a prefactor l/Air in which case e and g have to be replaced by \/Aite 
and y/Airg and the condition reads 2eg = n. 
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is locally coupled to the magnetic current j^x). The world line L of the electric charges 
are now represented by a gauge field 

P;„ = e^ x J XK (x;S) (193) 

with some irrelevant surface S whose boundary is L. This field energy is invariant under 
the dual charge gauge transformations 

K - K + K- ( 194 ) 

At higher temperatures, the monopole world lines proliferate on account of their entropy. 
In this phase, they can be described by a disorder field theory and the total action reads 

A = j d'x ji (F, v - F; v ) 2 + (d,0 - gA,)*} . (195) 

By integrating out the angular variables 6 of disorder, we are left with a transverse 
mass term of the A, field 

d'^Al 2 . (196) 

As a consequence of this, the previously irrelevant surfaces S of the charge gauge 
fields become massive and give rise to a continuing force between the boundary lines. 
Going back to the electromagnetic description, the action can be written as 

A = j d 4 x ji {F, v - f, v f + ij,A, - i A, {^-e, vXK d v f XK - j^j + Tj^-rJj} • ( 197 ) 

Here A, plays the role of a Lagrangian multiplyer which enforces j ^ to be a gauge field 
for the monopole world lines L. 

In the high-temperature phase, the monopoles proliferate and can be treated 
as an ordinary fluctuating transverse field giving rise to the mass term ^-A^. After 
this, Ap can be integrated out and produces a term l/2m 2 A (Ic^k^/sk) . The /^j, -field 
is responsible for the generation of the mass term of the surfaces S in the gauge field 
description of the electric charges, as can be verified by going again through a duality 
transformation. 

Near the critical temperature, the fluctuating j-currents can be represented by a 
disorder field theory of the |?/>| 4 -type and the total action reads 



A = j d 4 x \- (F^ - f^) + ijf.A, - iAi -t llv \ K d v f\ K - j f 



+i(^-zi,)vr+m 2 H 2 +!iv>i 4 }. 



198) 



ACKNOWLEDGMENTS 

I thank Dr. A. Schakel and M. Kiometzis for useful discussions. 



29 



[1] H. Kleinert, Gauge fields in Condensed Matter, 

Vol. I, Superflow and Vortex Lines, Disorder Fields, Phase Transitions , 

Vol. II, Stresses and Defects, Differential Geometry, Crystal Defects , World Scientific, 

Singapore, 1989. 

[2] H. Kleinert, Phys. Lett. A 130 (1988). 

[3] H. Kleinert and W. Miller, Phys. Rev. Lett. 56, 11 (1986); Phys. Rev. D38, 1239(1988). 

[4] Kleinert, H., Gauge Theory of Dislocations in Solids and Melting as a Meissner-Higgs 
Effect, Lett. Nuovo Cimento 35, 405 (1982) 

[5] S. Elitzur, Phys. Rev. D 12, 3978 (1975). 

[6] P. A.M. Dirac, Principles of Quantum Mechanics , 4th ed., Clarendon, Cambridge 1981, 
Section 80; 

[7] J. Schwinger, Phys. Rev. 115, 721 (1959); 127, 324 (1962). 

[8] M. Kiometzis, H. Kleinert, and A. Schakel, Phys. Rev. Lett, (in press) 

[9] M. Kiometzis and A. Schakel, Int. J. Mod. Phys. B 7, 4271 (1993). 

[10] H. Kleinert, Lett. Nuovo Cimento 35, 405 (1982). See also the more detailed discussion 
in Ref. [1], Vil I, Part 2, Section 13. The action (106) was derived there [see, in 
particular, Eq. (13.30)]. 

[11] M. Gabay and A. Kapitulnik, Phys. Rev. Lett. 71, 2138 (1993); 

[12] S.-C. Zhang, ibid., 2142 (1993); 

[13] M.T. Chen, J.M. Roessler, and J.M. Mochel, J. Low Temp. Phys. 89, 125 (1992). 

[14] D. Nelson, Phys. Rev. B 18, 2318 (1978); D. Nelson and B.I. Halperin, Phys. Rev. B 
19, 2457 (1979); A. P. Young, ibid., 1855 (1979). 

[15] D. Nelson, Phys. Rev. B 26, 269 (1982). 

[16] W. Janke and H. Kleinert, Phys. Lett. A 105, 134 (1984); Phys. Lett. A 114, 255 (1986); 
[17] W. Janke and H. Kleinert, Phys. Rev. Lett. 61, 2344 (1988). 



30 



